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Thesis of the author [T] stated that curvelet transform has non-linear approximation 
rate M~ 2 for functions of two variable that are C 5 apart from C 3 edges. This means 
that extra smoothness allowed to remove log-factor from well known approximation 
rate M~ 2 (log(M)) 3 that hold if function is C 2 apart from C 2 edges. Here the some 
results from thesis are generalized by replacing C 5 smoothness assumption by C 3 
assumption. For definitions, look pp. 

Theorem 1 Let f e F^ yn , N > n, p e S be point that minimizes L = \D a Rg(b — p)\ 
and 6' ~ ka 1 ^ 2 be the angle between major axis of 'jaQb and tangent of S at point p. 
Then for any K > and < e < 2 



Proof. First we assume that L < e l 2 . With this assumption the point p is on 
major axis of 7 (see proof of theorem 15 in pQ). 

Let assume that for regions Ri C 1R 2 holds 




1*1 < « 1/2 

n > 3, |^| > a 1 / 2 , Ik] 1 -^ 2 > L 
n>3, \9'\ > a 1 / 2 , l^ 1 ^ 2 < L 



(1) 



U 3 i= _ 2 Ri = R 2 and 



(2) 



and define 



fi 



= IrJ. 



(3) 



Then 





Next we assume that for regions Ri C 1R 2 hold 



Rif\Ri = 0. 



(5) 



/ 

JR- 



(6) 



With any functions P{ we can write 

laObfi 

jr 2 

= / Ja8b(lRif ~ ^-RiPi) + / Jaeb^RiPi 
JR 2 JR. 2 

= [ 7«»(W-Vi)+/ lae b {l R A + ^RPi-^RPi) 

JR. 2 JR. 2 
= / la8b(lRj - IrA) + / lae b l Rt yjRPi ~ I 7^1^) 

Jr. 2 Jr. 2 Jr. 2 

= / la8b{f ~ Pi)+ _ 7a0bP - / 7a0b-Pi 

Now, for i = —2, . . . , 2, we define regions Ri and i?j as illustrated in Figure [l] and 

P 3 = \U 2 = _ 2 ^. (7) 



Regions i? 0; -Ri an d R2 are also illustrated with greater details in Figures |2j [3j |4j On 
those figures lengths d and ft, are as follows; 

d « a 1/2 £r\ /i « afc" £ . (8) 

Also is assumed to be small, i.e. sin(6 |/ ) ~ 6'. Case of "big" angles is omitted here 
since acceptable estimates are very straightforward to produce in that case (one can 
for example take d ~ a and reproduce the rest of proof quite straight). 

First, in the case i = 3, by rapid decay of 7, we get 

a 3/4 

We concentrate next to case i = 2. If we define P 2 such that it is polynomial in 
direction of minor axis of "f a e b , then, since vanishing moments of 7 a efe, 

/ lcmP2 = 0. (10) 
JR2UR2 

Moreover, because rapid decay of 7 and dimensions of R2 (look figure [2]), 

a 3 / 4 



laQ b P2 
R2 



< ~ (11) 



By choosing P 2 so that slices of P2 in the direction of minor axis of 7^ are always 
second order Taylor polynomials, developed at major axis of j a gb for slices of /, we 
get 



laBbf 

R 2 



< a 3 / 4+3 . (12) 




Figure 1: Illustration of regions Ri. 
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Figure 2: (a) Illustration of regions R2 and i?2- 




Figure 3: Illustration of regions R\ and R\. 
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Figure 4: Illustration of regions R and R . 



Now we look the case i — 1. In this case we again define Pi by using second order 
Taylor polynomials. However this time we define those along radial lines that intersect 
in point Q in Figure |3j Points where polynomials are developed can be at major axis 

Of 7a0&- 

By doing integration in the polar coordinates r and a (origo at point Q) we get 

laebPi = 0. (13) 



/ . 



This is because when when integrating first respect to r we get zero (because of 
vanishing moments of jadb- 



By looking shape of Ri and using rapid decay of 7^ we get again 



ladbPl 



Ri 



< 



,3/4 



l+ \ k Ul-e/2)K L K- 



(14) 



When making estimate for 



S Rl laebU - Pi) 



we divide R\ into subregions R\i like 



illustrated in figure XX. First we notice that because of the way how we defined Pi, 

xeR^ => \f(x)-P 1 (x)\<(ld) 3 . (15) 
Also area m(Ri j i) of region Rij is bounded by 

m{R u ) < m{R m i- s/2 ) « as 2 w a^ 2 k~ £ ' 2 . (16) 
By using these two estimate and rapid decay of 7, we get 



la9b{f - Pi] 



Ri 



ll-e/2 



< 



E 



-3/4 



i,« l + \l\ K L* 



< 



Z=-|fc| 1 -/ 2 
a 3/2 k -e/2 d 3 



{Idf 



(17) 



1 + L K 

a 3/2-3/4+3/2^,-3-E/2 

1 + L K 

This is clearly acceptable. 

Cases i = — 1 and i = —2 are identical to cases i = 1 and % = 2, so those are omitted. 

Now only the case % = is left. But Rq is exactly the same region as Rt in [1] and it 
can be handled similar techniques as in there. For the sake of completeness we give 
here the major steps (and also correct some minor missprints). 

The major point is that Rq is heavily "corrupted" by discontinuity curve S, i.e. we 
cant apply Taylor polynomials directly. Instead of that, we do such a change of 
variable that "twist" S to straight line inside region Rq. 
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First we define x\- and 22-coordinate axes as illustrated in Figure |4} Inside Ro, with 
small enough scales a, curve S can be considered as function g{x\). 



We define the twisting operator T g : R — >■ Ro by formula 

Tg(X 1 ,X 2 ) : = 

Notice that 



(x u x 2 + >L -j^-g(x 1 )) 
(x 1 ,x 2 ) 



xi,x 2 ) e -Ro 
x 1 ,x 2 ) i Ro 



;is) 



(19) 



T g Ro — Ro 

and if apart from S, f is three times continuously differentiable (with bounded deriva- 
tives) and first three derivatives of g are continuous and bounded, then first 3 deriva- 
tives of the function 

f(x) := f(Tgx), (20) 
in direction of Xi-axis are bounded and continuous inside Ro. 
Change of variable 

x = T g y, (21) 



gives 



where 



and 



Since 



dx 



det( J(y))dy , y e R 
dy , y £ Ro 

^ 2 dy , yeRo 
dy , y £ R 



(1- 
dy 



sqn(y 2 )^fp-)dy 



yeR 
y^Ro 



f{x)j a b (x)dx 



Ro 



f(y) det(J(y))% eb {y)dy 



Ro 



h{y) :=/(y)det(J(y)) 



f(y)det(J(y)) , y e R 



f(y) 

laeb{y) ■= ^ a eb{T g y). 

g{yi)\^yli \g'{yi)\ ^ M, °< 



y i Ro 



h - \y 2 \ 
h 



< 1, 



it is quite clear that 



l%)l<i, ye 



dh{y) 



dyi 



< 



1 2/1 1 „ cl R 

~h ■, y e n 



y^Ro 



(22) 



(23) 

(24) 
(25) 

(26) 
(27) 
(28) 



S 



and for 2 < m < 3 



d m h(y) 



dyf 



1 , y ^ -Ro 



(29) 



On the border of Rq the function %e b is discontinuous but all decay properties of %e b 
(and it's derivatives) remain. Also 



i(y) = i(y),Vy Rq 



(30) 



However, unlike j a ob, the function %g b does not have directional vanishing moments. 
Because of that we will "recreate" function 7^: 



f(x)~f a g b (x)dx 

Ro 

= / h(y)% eb (y)dy 

JRo 

= I Hy)(%0b(y) - la6b(y) + laeb(y))dy 
JRo 



(31) 



Ro 



h(y)(%0b(y) - la8b{y))dy + / h-f a9b (y)dy. 

JRo 



The last of these two integral is handled now very similarly what we did with re- 
gion R 2 : we create polynomial P by defining it as second order Taylor polynomials 
separately for each slice of function h that are aligned to rrx-axis. This way we get 



/ h-f a9b = / -f a g b (h -P )+ / laebPo ~ I laebPo- 

J Rq J Rn JRaURo JRo 



(32) 



'-Ro JR0UR0 
Because of vanishing moments of 7 we have 



/ 

JR( 



laObPo = 0, 



(33) 



iRoURo 

and because of rapid decay of 7 we have 

< 



/ 7a 


9ft -Pq 


JRo 





,3/4 



1 + \k\^+ £ / 2 ) K L K ' 



(34) 



By using regularity of function h, rapid decay of 7 and dimensions of R we get with 
straightforward calculation that 



/ ladb{h - P ) 
JRo 



rs ~ / I , I 7 1 4 



-3/4 



1 + WL K 



(35) 



and by remembering also that 7 is differentiable (infinitely many times), we get 

a 3 / 4 



/ 

JRo 



Ky){iaeb{y) -iaBbiy))dy 



< 



1 + \ k \ 3+£ l k 



(36) 



Details for these calculations can be founded from proof of theorem 14 in [Tj (only 
difference is that here we already write factor L K , that follows from rapid decay of 
7, visible already here). 

Finally we turn to case L > /c 1_£ / 2 . Here we get straightforwardly estimate 



1a9bf 
Ro 



a 3 / 4+3 a 3 / 4 

~ - max ^ i +fc (ww a } (37) 



by using similar techniques as with case i = 2 before. Notice that all exponents related 
to K can be simplified to form used in theorem since there is not any limitation for 
K and 1 — e/2 > 0. ■ 

Theorem 2 Let fuc be M-term non-linear approximation of f by using curvelets. 
Iff£F 3 , 3 ,then\\f-U w \\ 2 2 <0(M- 2 ). 

Proof. Proof is exactly same as proof of Theorem 17 in pQ, now we just can apply 
the improved decay estimate ([T]). ■ 
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